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Amperes law states that the magnetic moment of a ring is given by current times the area enclosed. 
Also from equilibrium statistical mechanics it is known that magnetic moment is the derivative of 
free energy with respect to magnetic field. In this work we analyze a quantum double ring system 
interacting with a reservoir. A simple S-Matrix model is used for system-reservoir coupling. We see 
complete agreement between the aforesaid two definitions when coupling between system and reser- 
voir is weak, increasing the strength of coupling parameter however leads to disagreement between 
the two. Thereby signifying the important role played by the coupling parameter in mesoscopic 
systems. 
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I. INTRODUCTION 

Mesoscopics, a fertile branch of Physics deals with sys- 
tems and phenomena in the scale of nano- to microme- 
ters. The distinguishing feature of these systems is that 
at extremely low temperatures an electron retains it's 
phase coherence throughout the sample. These systems 
have revealed a new range of unexpected quantum phe- 
nomena, often counter-intuitivei*^. The notion of the 
usual ensemble averaged transport coefficients such as 
the conductivity, i.e., local and material specific, has 
to be replaced by that of conductance, i.e., global and 
operationally specific to the sample as well as the na- 
ture of probes of measurement. Some novel and hitherto 
unheard of features in classical physics, e.g., non-local 
current- voltage relationships 3 , breakdown of Ohm's law, 
absolute negative resistance (four probe)^, normal-state 
Aharonov-Bohm effect, quantization of point-contact 
conductance (Landauer formalism) 5 , universal conduc- 
tance fluctuations (new form of ergodicity), persistent 
currents 6 , spin-polarized transport-' 8 , coulomb blockade 9 
and many novel effects arising due to electron correla- 
tions, have been observed in these systems. Interpreta- 
tion of these require full recognition of the wave nature of 
quasi particles and keeping track of their phase coherence 
over the entire sample including the measurement leads 
and probes (quantum measurement process). Even the 
equilibrium properties, are very sensitive to the nature of 
statistical ensemble used. The results differ qualitatively 
from one ensemble to anotheri 



II. MOTIVATION 

Recently it has been proposed that these systems will 
provide a testing ground for verifying the violation in 
the basic laws of thermodynamics^. This behavior has 
been explained by taking recourse to the effects of en- 
tanglement, through which the quantum system is so in- 
terlinked with the bath that the resulting behavior of 



the system alone cannot be treated within a conven- 
tional thermodynamic approach. Here the finite cou- 
pling between the bath and the system plays a crucial 
role. It should be noted that equilibrium thermody- 
namics of the super-system comprising of system (sub- 
system) plus bath, does not imply standard equilibrium 
thermodynamics for the sub-system alone. In-fact the 
thermodynamic equilibrium properties of the system de- 
pend on the coupling parameter, unlike equilibrium sta- 
tistical mechanics. For example, it has been shown via an 
analytical treatment that the mean orbital magnetic mo- 
ment, a thermodynamic property, is determined by the 
electrical resistivity (which is related to system-bath cou- 
pling parameter) of the materialii. What is crucial for 
dissipative diamagnetism is that system-bath interaction 
has to be treated exactly, there is no clear cut separa- 
tion between what is the system and what is the bath- 
both are inexorably linked to one many-body system. In 
this work, motivated by the above results, we provide a 
simple example wherein the equilibrium properties are 
determined by the system-reservoir coupling parameter 
in a non-trivial manner. Our results follow from the 
consideration of the dephasing of a single particle quan- 
tum coherence while the earlier dynamical treatments re- 
quired in addition to quantum coherence, entanglement 
between system and bath. It is in this spirit that we 
study the persistent current densities in a quantum dou- 
ble ring system coupled to a reservoir via a simple volt- 
age probe method due to Biittiker— . We explicitly show 
that when the coupling parameter is very small there is 
perfect agreement between the magnetic moments calcu- 
lated from the local currents (via Amperes law) and that 
from the derivative of free energy with respect to mag- 
netic field, increasing the strength of coupling parameter 
however leads to disagreement between the two. 



III. BACKGROUND 

It is well known that spontaneous currents which never 
decay can flow in super-conducting systems. In 1983, 
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Biittiker , Imry and Landauer— first predicted that, a 
normal metal ring threaded by an Aharonov-Bohm flux 
cj>, in the phase coherent domain carries persistent cur- 
rents . These arise because the magnetic flux breaks the 
time reversal symmetry thus inducing currents. This is 
a quantum effect and the total current flowing in the 
ring is related to the derivative of free energy with re- 
spect to flujii&ii. For clean rings it is periodic in flux 
with period of 27T0/0O- The magnetic flux incidentally 
plays the same role as a periodic potential in the Bloch 
sense, and therefore one gets the band structure, elas- 
tic scattering if present in the loop induces gaps in the 
spectrum and reduces the magnitude of persistent cur- 
rents. Later in 1985, Biittikei*^ investigated the effect of 
a reservoir coupled to a ring. This simply means break- 
ing the phase coherence. Electrons enter the reservoir 
lose their phase memory and are re-injected again from 
the reservoir with an uncorrelated phase. The S-Matrix 
for coupling between system and reservoir is given by- 

-(a + b) Ve \/e 
S., = | a b 

y/e b a 



IV. MODEL 

Let us consider the double ring system as shown in 
FIG. 1. The static localized flux piercing the loop is nec- 
essary to break the time reversal symmetry and induce a 
persistent current in the system. The geometry we con- 
sider is a one-dimensional ring with an attached bubble 
and a lead connected to a reservoir at chemical potential 
/j. The reservoir acts as an inelastic scatterer and as a 
source of energy dissipation. All the scattering processes 
in the leads including the loop are assumed to be elas- 
tic. Hence there is complete spatial separation between 
the elastic and inelastic processes. The loops JlJ2aJ3Jl 
and JlJ2bJ3Jl enclose the localized flux $. However, the 
bubble J2aJ3bJ2 does not enclose the flux <!>. The same 
S-Matrix coupler couples the double ring to the reser- 
voir, i.e., = Sj, for the other two junctions we take 
symmetric couplers 2 ^. 



Sj2 = Sj3 



wherein, a = ±(V1 - 2e - 1) and b = ±(V1 - 2e+ 1), 
for e — * the system and reservoir are decoupled while 
for e — > 0.5 the system and reservoir are strongly cou- 
pled. The above S-Matrix satisfies the conservation of 
current™, and accounts for the possibility of strong to 
weakly coupled reservoirs through the coupler "e" . One 
of the important conclusions of the work was that the 
magnitude of persistent current flowing in the loop de- 
creases with increasing coupling strength e but without 
any change in their nature, i.e., diamagnetic or param- 
agnetic. This effect is solely due to exchange of carriers 
between reservoir and ring (dephasing), and also this is 
true if the lead connecting reservoir to ring has a charging 
energy much less than the level spacing. Experimentally 
persistent currents in both open and closed systems have 
been observed^iSiiLiS and these observations have given 
rise to a spurt in theoretical activities^. 




FIG. 1: One dimensional mesoscopic ring coupled to a bubble 
with a lead connected to a reservoir at chemical potential /i. 
The localized flux <j) penetrates the ring. 



The waves in the four arms of the system depicted in 
FIG. 1 are related as follows: The waves incident into 
the branches of the double ring system are related by the 
S-Matrices for J1J2 arm by- 



for J2bJ3 arm by- 
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Here kli, M2, kl% and fc/4 are the phase increments of 
the wave function in the absence of flux, 01, 0:2,0:3, and 
0:4 are the phase shifts due to flux in the arms of the con- 
sidered double ring system. Clearly, 01 + 02 + 04 = 
where $ is the flux piercing the loop and $0 is the flux 
quantum also 01+03 + 04 = thus 02 = oj 21 ' 22 . 
The current densities (in dimensionless form) in the 
various arms of the system are given as follow o 12 ! 21 

h = |zi| 2 - M 2 , h = \X2? - |y 2 | 2 , h = M 2 - bl 2 , 

I4 = — |rc4 1 2 + |?/4| 2 , wherein complex amplitude for prop- 
agating waves x\, X2,X3, X4, 2/1,2/2; 2/3, 2/4, are as depicted 
in Fig. 1. The induced current densities in the various 
arms of the loop are assigned labels 1\ , I2, 13 and I4, while 
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the arm lengths of various parts of the ring are h,h, h 
and I4. 

Now a question may arise as to why we consider 
the persistent current densities in this special geometry 
of a double ring. This geometry exhibits the current 
magnification effect in equilibrium^ which is a quan- 
tum phenomena. The current magnification effect, a 
purely quantum phenomena arises in equilibrium and 
non-equilibrium systems^S&Si and leads to large orbital 
magnetic moments. A double ring system 2 ^, which does 
not exhibit current magnification effect will show no dis- 
agreement between the magnetic moments calculated 
from the local current densities (via Amperes law) and 
that from the derivative of the free energy, which we have 
separately verified^. 



V. THEORY 

We know that magnetic moment is the derivative of 
free energy with respect to magnetic field, /i = — 
wherein F is the free energy and H is the magnetic field 
enclosed by the system. From text books we know that a 
current carrying loop behaves as a magnet, in other words 
Amperes law which states- magnetic moment is current 
multiplied by the area enclosed. The orbital magnetic 
moment density for a system coupled to reservoir can 
be calculated in two ways, first from the formulation of 
Akkermans, et alj2L2& one can calculate it as follows- 

1 djlndetS] 
d V = 7T- TT7 dE (1 

2m o<p 

Here, is the differential contribution to the mag- 
netic moment at energy E, and S is the on-shell scat- 
tering matrix. In the system considered in Fig. 1, the 
on-shell scattering matrix is just the complex reflection 
amplitude r. Also from Amperes law one can calculate 
the orbital magnetic moment density via the local cur- 
rents in the system. The orbital magnetic moment den- 
sity defined via local currents in a loop, depends strongly 
on the topology of the system, whereas the magnetic mo- 
ment densities calculated from the eigen spectrum as also 
in the formulation of Akkermans, et al, do not. This 
is special to only one-dimensional geometry where the 
eigen-spectrum is independent of variation in topology 
of the system . In fact there are infinitely many topo- 
logical structures possibl e) 22 ' 2 ? . If we consider our sys- 
tem as depicted in Fig. 1, to be planar and lying in 
the x-y plane then the magnetic moment density (/zi) 
can be viewed as being generated by current density I\ 
enclosing an area A r and by current density ^3 enclos- 
ing area Af,, i.e., [i\ — ^{I\A r + I^Ai,), wherein A r and 
Ai, are the areas enclosed by the ring (JlJ2aJ3Jl) and 
the bubble {J 2a, Jib, J 2) respectively. Another orienta- 
tion of the system in which the arm J2bJ3 is in the x-z 
plane gives /12 = \(I\A r + l2A r )/2 — I c A r /2 wherein 



I c = Ii + I2 is said to be the most appropriate general- 
ization of the equilibrium persistent current and which is 
consistent with Eq. 1, see [2j| for further details. Sev- 
eral other orientations are possible, for example, if the 
bubble lies in x-y plane and the ring lies in x-z plane, 
then [i z = ^(I^Ar — l2A r )/2 and [i y — ^l\A r . Even 
when our system lies in the x-y plane for fixed Zi , , ^3 5 
by deforming their shapes we can have different values 
of magnetic moment density along the z direction. It is 
also worth mentioning that the total magnetic moment 
(at temperature T — 0) of a representative system is 
obtained by integrating the magnetic moment densities 
up-to the Fermi wave- vector kf. 



VI. RESULTS 

From the Amperes law we calculate the orbital mag- 
netic moment density for some length parameters. After 
calculating we plot the orbital magnetic moment densi- 
ties obtained for the case of ^2 = -(hA r + l2A r )/2 = 
I c A r /2. In Fig. 2(a) we depict the plot of the energy- 
eigen values (normalised by it 2 ) of the closed system, and 
in Fig. 2(b) we plot the dimensionless orbital magnetic 
moment density fi2 obtained via the local persistent cur- 
rent densities as a function of the dimensionless Fermi 
wave-vector kl for different coupling parameters. Now 
lets test the equivalence of these definitions as e — > 0, 
i.e., the reservoir and system are almost decoupled. We 
see complete agreement. It can be noted from Fig. 2(a) 
that the ground state carries diamagnetic current while 
1st, 2nd and 3rd excited states carry paramagnetic cur- 
rent for small values of flux (which is obvious from their 
slopes), the 4th excited state carries diamagnetic currents 
while the fifth paramagnetic current. Similarly from the 
top most panel of Fig. 2(b), the e — > limit, it can be seen 
that the ground state carries diamagnetic current (mag- 
nitude is negative) while 1st, 2nd and 3rd excited states 
carry paramagnetic currents (magnitude is positive), the 
4th carries diamagnetic current while the fifth paramag- 
netic current. Indeed the two definitions are completely 
equivalent as far as the e — > limit (henceforth referred 
to as the weak coupling case) is considered. Now as we 
increase e, we notice a dramatic change from the weak 
coupling case (see middle panel in FIG. 2(b)). There 
are paramagnetic-diamagnetic jumps at those Fermi en- 
ergies wherein one would have expected pure diamagnetic 
or paramagnetic current. This behavior is more seen as 
one increases the coupling till the maximum (e = 0.5) is 
reached, see lowest panel of Fig. 2(b). Looking closely at 
the middle panel of Fig. 2(b), one can notice that the 
peaks broaden, noticeably the ground, 1st, 3rd and 5th, 
while the nature of the currents carried at 2nd and 4th 
levels changes qualitatively. Also, as one approaches 
the e = 0.5 limit the 1st, 3rd and 5th levels disappear 
completely, these are some of the other qualitative fea- 
tures which distinguish the cases depicted in Fig. 2(a) 
and (b) . This reaffirms that predictions from equilibrium 
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FIG. 2: Plot of (a) energy levels and (b) orbital magnetic moment density /j.2, for length parameters h/l = Ia/1 = 0.375, 
fo/l = 0.15, I3/I — 0.85 and flux = 0.1. In (b) the upper, middle and lower panels are for coupling strengths e — > 0, e = 0.05, 
and e = 0.5 



thermodynamics are in general not valid for a quantum 
system strongly coupled to a bath. In the conventional 
equilibrium treatment the properties of the bath appear 
only through a single parameter namely, the tempera- 
ture T and nowhere the coupling parameter appears in 
the magnitude of equilibrium physical quantities. Our re- 
sults reaffirm that equilibrium properties are determined 
by the coupling parameter in the quantum domain. The 
finite coupling can lead to qualitative changes (and not 
just the broadening of levels or persistent current peaks) 
from that of the predictions of equilibrium statistical me- 
chanics as shown above. 

We have verified separately that the orbital magnetic 
moment density calculated by the formulation of Akker- 
mans, et al, (Eq. 1) is similar to that in Fig. 2(b), for the 
geometry considered hereSi. In addition we also have 
calculated the orbital magnetic moment density for dif- 
ferent topological configurations, e.g., (Xi, (x z , fi y (see the 
Theory section for further details) and obtained results 
not in consonance with that calculated from the eigen- 



spcctrum (equilibrium statistical mechanics), thus bol- 
stering the fact that the orbital magnetic moment den- 
sity calculated from the local current densities is inher- 
ently linked to the topology of the system. As already 
pointed out, the orbital magnetic moment density calcu- 
lated through the local current densities is qualitatively 
different (nature) from that calculated from the closed 
system energy eigen-spectrum which is independent of 
topological variations. 

We have seen these novel features not only for one 
set of length parameters but for many different set of 
length parameters. In Fig. 3(a), we plot the energy 
eigenvalues (normalised by 7r 2 ) as a function of flux, and 
in Fig. 3(b) the dimensionlcss orbital magnetic moment 
densities as a function of the dimensionlcss Fermi wave- 
vector kl for different length parameters. Herein also 
it can be seen from Fig. 3(a) that the ground and the 
first excited states carry diamagnetic currents while 2nd 
and 3rd carry paramagnetic currents for small values of 
flux (which is obvious from their slopes), the 4th and 
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FIG. 3: Plot of (a) energy levels and (b) orbital magnetic moment density /j.2, for length parameters h/l = Ia/1 = 0.375, 
I2/I = 0.05, I3/I — 0.95 and flux = 0.1. In (b) the upper, middle and lower panels are for coupling strengths e — > 0, e = 0.05, 
and e = 0.5 



5th again carry diamagnetic currents. Similarly from the 
top most panel of Fig. 3(b), the e — > limit, it can be 
seen that the ground and first excited state carry dia- 
magnetic current (magnitude is negative) while 2nd and 
3rd carry paramagnetic currents (magnitude is positive), 
the 4th and 5th carry diamagnetic currents. The two 
definitions are again completely equivalent as far as the 
e — > limit is considered. Again as we increase e, we 
notice the same change from the weak coupling case (see 
middle panel in FIG. 3(b)). There are paramagnetic- 
diamagnetic jumps at those Fermi energies wherein one 
would have expected pure diamagnetic or paramagnetic 
current. This behavior is more seen as one increases the 
coupling till the maximum (e = 0.5) is reached, see lowest 
panel of Fig. 3(b). Thus, in contrast, to the case of a sin- 
gle quantum ring coupled to a reservoir as was considered 
by Biittiker, wherein coupling only led to a broadening 
of energy levels^, in the case of a quantum double ring 
system considered here, in addition to level broadening 
one also sees a change in nature of currents as one 



increases the strength of coupling to the reservoir. 

One should also mention here that there is a drawback 
in modeling the inelastic effects by this model, if one 
changes the position of attachment of lead and the dou- 
ble ring system. Qualitatively, different results for the 
orbital magnetic moment density would be obtained as 
by definition it involves the length parameters of our sys- 
tem, These length parameters will of-course change with 
position of junction (Jl) . Only in the very weakly cou- 
pled regime can the specific position of lead to reservoir 
attachment be ignored. Of-course this drawback does not 
exist for the system investigated by Biittiker in Refill 
wherein the orbital magnetic moment density would be 
same irrespective of the position where lead is attached. 
However, this simple model for coupling between system 
and reservoir is enough in order to bring out the impor- 
tance of finite coupling between system and bath, vis a 
vis equilibrium thermodynamics. 
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VII. CONCLUSIONS 

In this work we have shown using a very simple model 
of system-reservoir coupling that the equilibrium proper- 
ties are determined by the strength of coupling parame- 
ter. This is consistent with the recent findingsi2*ii, that 
in the quantum domain, equilibrium properties of the 
sub-system are related to dissipative coefficients arising 
due to subsystem-bath coupling, unlike the predictions of 
conventional statistical mechanics. These fully dynami- 
cal studies in the quantum domain invoke coherence and 
entanglement between system and bath degrees of free- 
dom to come to the same conclusion. However our treat- 
ment is simple and invokes the single particle coherence 
which is disrupted by the presence of the reservoir (with- 
out bringing the notion of quantum entanglement be- 
tween system and bath). However these finite coupling 
induced qualitative changes can be observed only in hy- 
brid rings which exhibit current magnification effect, a 



purely quantum phenomena, in equilibrium 21 . We have 
verified separately that multiple ring structures which 
do not exhibit current magnification effect do not show 
these qualitative changes apart from features expected 
from broadening of energy levels^. These results can be 
verified experimentally by attaching a voltage probe cou- 
pled to a nanoscopic semi-conducting system as depicted 
in Fig. 1. Here the voltage probe serves the purpose of a 
reservoir. One can change the parameters of the junction 
between system and lead by applying appropriate gate 
voltage (to simulate the effects of coupling parameter) 
beneath the junction and measure the orbital magnetic 
response as function of the gate voltage. 
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